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ABSTRACT 

The morphology of the atomic hydrogen (H i) disk of a spiral galaxy is the first com- 
ponent to be disturbed by a gravitational interaction such as a merger between two 
galaxies. We use a simple parametrisation of the morphology of H i column density 
maps of Westerbork Hi Spiral Project (WHISP) to select those galaxies that are 
likely undergoing a significant interaction. Merging galaxies occupy a particular part 
of parameter space defined by Asymmetry (A), the relative contribution of the 20% 
brightest pixels to the second order moment of the column density map {M20) and 
the distribution of the second order moment over all the pixels (Gm)- 

Based on their Hi morphology, we find that 13 % of the WHISP galaxies are in 
an interaction (Concentration- M20) and only 7% based on close companions in the 
data-cube. This apparent discrepancy can be attributed to the difference in visibility 
time scales: mergers are identifiable as close pairs for 0.5 Gyr but ~ 1 Gyr by their 
disturbed H I morphology. Expressed as volume merger rates, the two estimates agree 
very well: 7 and 6.8 x 10"'^ mergers Gyr~^ Mpc~'^ for paired and morphologically 
disturbed H i disks respectively. 

The consistency of our merger fractions to those published for bigger surveys such 
as the Sloan Digital Sky Survey, shows that H i morphology can be a very viable way 
to identify mergers in large H i survey. The relatively high value for the volume merger 
rate may be a bias in the selection or WHISP volume. The expected boon in high- 
resolution Hi data by the planned MeerKAT, ASKAP and WSRT/APERTIF radio 
observatories will reveal the importance of mergers in the local Universe and, with the 
advent of SKA, over cosmic times. 
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1 INTRODUCTION 

Mergers of galaxies is a driving factor in their evolution over 
cosmic times. Several schemes to identify merging galaxy 
pairs have been developed in the past decade, many based 
on the number of physically close pairs (in both sky co- 
ordinates as well as redshift) or on the characterization of 
the disturbed appearance of galaxies due to gravitational 
interaction (often through visual i n spection, e.g.. lArDlll973! : 
IVorontsov-Velvaminov et all I2001I : iDarg et al.1 12009| ). Both 
these techniques have been used to determine the interac- 
tion fraction in the local universe as well as out to high red- 
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shift in Hubble images. Using N-body simulations, one can 
determine how long a merger will be identified as such by 
both techniques; the galaxy pair is close enough, the galaxies 
look sufficiently disturbed. 

Because the volume probed increases with redshift, 
there were until recently -paradoxically- better measures of 
the interaction fraction for higher redshift than for the local 
Universe. The S PSS search for clos e (and disturbed looking) 
pairs of galaxies (|Darg et al.ll2009l l added the valuable local 
Universe inte r action fraction, improving on the estimate by 
iPatton et al.l |l997 ), from close pairs in the Uppsala Gen- 
eral Catalog (|Nilson|[l973l ). A compilation of merger frac- 
tions determined as a function of redshift is shown in Figure 
[T]and Table [1] The scatter in the merger fractions, even for 
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Table 1. The reference, data-set, and method for merger fractions 
in the local and distant Universe. 

Reference Data-set Criteria 

Morphology 



Conselice et al. (20031 
Conselice et al. (20081 
Conselice et al. (20091 
Conselice et al. (2005 1 
Lotz et al. (2008J) 
Scarlata et al. (20071 



HDF-N CAS 

HUDF CAS 

ECS, COSMOS CAS 

HDF-S CAS 

ECS G/K'ho 

COSMOS CAS-fG/Af20 



Pair statistics: 



Lin et al. (20041 
Lin et al. (20081 
Kartaltepe et al. (2007) 
de Ravel et al. (2009 1 
^assate^t_^l^ (2005) 
Le Fevre et al. (20001 
Patton et al. (19971 
Patton et al. (20021 
De Propris et al. (2007a1 



DEEP2 

DEEP2 

COSMOS 

VLT/DEEP 

GOODS 

CFRS, HST 

CNOCl 

CNOC2 

MGC 



"T ' — 

o Conse 
• Conse 
o Conse 
■ Conse 
□ Conse 



ice et a 
ice et a 
ice et a 
ice et a 
ice et a 



(2003b) 
2003b 
2008b 
2009b 

2006) HDF-S 



* Lotz et al. , 

' Scarlata efai. (2007) 

® Darg et al. (2009) 



o Lin et al. (2005) 

• Kartaltepe et al. (2007) 
o Lin et a! (2008) 

• Lin et al. 2005) 

n Cassata et al. (^2005) 
■ Le Fevre et al. (2000) 
de Ravel et al. 2008) 
' de Propris et al. (2007) 
' Patton etal. (1997) 
V Patton et al. 2002 



those determined over the same data is striking and it is 
tied to the definition of what constitutes a pair or disturbed 
morp hology (see for instance the discussion in iGenel et al.l 
I2OO9I I. The definition of morphologically disturbed became 
quantified in several schemes of morphologic al param eteri- 
sation schemes (e.g.. IConselicell2003l : iLotz et alfl2004 l. Ob- 
servational uncertainties are the time scale a merger is iden- 
tifiable as one, the completeness of the various samples for 
each technique and the volumes considered. Similarly, a cur- 
rent substantial theoretical effort is to map the dark mat- 
ter halo merger rates onto actual observable galaxy mergers 
l|Hopkins et al.|[201(]| . and reference therein). 

The morphological studies are largely based on opti- 
cal, mostly B-band and restframe UV data. The reasoning 
goes that mergers trigger star-formation and the resulting 
increased surface brightness make the disturbed morphology 
easier to pick up (although the increa se in star-formation 
is not a given, see iRobaina et"alll2009l 1. However, with the 
emergence of new and refurbished radio observatories in 
preparation for the futu re Square Kilometre Array (SKA; 
ICariUi fc Rawlingsll2004l l. a new window on merger rates 
over cosmic times will be opening up: the 21 cm emission 
line of atomic hydrogen gas (Hi). The two SKA precur- 
sors, the South Af r ican Karoo Array Telescope (Mee rKAT; 
iBooth et al] l2009l : Ijonad l2007l : Ide Blok et al] |2009|) and 
the Austrahan SKA Pathfin de r (AS KAP; |johnstonlT2007l : 
IJohnston erahl |2007| . l2008al lbl. |2009| ) stand poised to ob- 
serve a large number of Southern Hemisphere galaxies in H I 
in the nearby Universe (z< 0.2). In addit ion, the Extended 
Very Large Array (EVLA; iNapieJ I2OO6I I and the APER- 
ture Tile In Focus instrume nt (APERTIF: IVerheiien et al] 
l2008l : lOosterloo et al.1 l2009l l on the Westerbork Synthesis 
Radio Telescope (WSRT) will do the same for the North- 
ern Hemisphere. The advantage of Hi observations is that 
it contains both morphological and kinematic information 
of spiral disks. There is ample anecdotal evidence of dis- 
turbed Hi morphology during a merger (see the compila- 



Figure 1. The merger fraction {fmerg) as a function of red- 
shift (z). The black points are based on quantified morphology 
estimates, the gray points are based on galaxy pair counts. The 
Conseli ce et al papers are based purely on the CAS classification 
system . iLotz et al ] (|200d l uses the Gini/M20 classification and 
IScarlata et al.l | |2007^ both the CAS a s the Gini and M20 param- 
eters. The results by iDe Propris_et_a l. { 2007a) an d [Darg: et al.l 
('2009') are hybrid approaches: iDe Propris et al.l (I2007ah l o oks at 
pair statistics but also galaxy asymmetry and lPar^Tt al.l l|2009l l 
used the vis ual identification of a merging pair in the Galaxy- 
ZOO proiect lLiiitott fc the GalaxyZOO teami l|2009l l. See for the 
data-sets used and the references Table [Tl 



tion in iHibbard etaLlbOOlf FI In this series of papers, we ex- 
plore primarily the signature of gravitational interaction on 
the morphology of the (face-on) Hi disk. This is a suitable 
complement to any kinematic signature, which will be most 
clear in edge-on disks. Our motivation to move to the H I 
perspective is that (a) the gas will be disturbed before the 
stellar disk is, (b) the H I morphology will be more sensitive 
to minor interactions, which may dominate the number of 
interactions and (c) the H I morphology will be intrinsically 
sensitive to gas-rich interactions. Minor and gas-rich inter- 
actions are expected to dominate at higher redshift, which 
makes an Hi perspective at low redshift a good local com- 
parison. 

In the previous papers in this series, we com- 
pared the Hi morph o logy t o that in other wavelengths 
(|Holwerda et all |2009| . l2011bh and found it to be at least 
as good a tracer of mergers as any other wavelength. We de- 
fined an Hi parameter space to identify interacting galaxies 
iHolwerda et aLll201ld l and derived a time-scale for inter- 
actions to reside in this param eter space (this paper's com- 
panion, |Hol^rd^^e£|al][201l3)- In this paper, the aim is to 
combine the morphological identification of mergers with the 
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Figure 2. bottom panel: The distribution of Hubble types in 
WHISP (black line) and the part of the UGC that conforms to the 
selection criteria for WHISP (thick gray line), top panel: per- 
centage of the UGC catalogue observed in WHISP as a function 
of type. There is a clear preference in the WHISP selection for 
later-type galaxies but Ellipticals are not specifically excluded and 
large early types would make the H I flu x cut. Hubble type deter - 
minations are from the 2MASS survey i|Ki einmann etal .111994) . 
not the UGC. 



time-scales into a merger fraction and rate for the WHISP 
sample. The organisation of this paper is as follows: in sec- 
tion 15.21 we briefly describe the morphological parameters 
and selection criteria, in section |3] we discuss the limitations 
and applicability of these in the context of H idata, in sec- 
tion UJ sections 14. II and 14.21 describe the WHISP basic data 
and H I column density maps. In section 14.31 we derive the 
volume representative for the WHISP survey. In section [5l 
we derive the merger fraction based on the number of pairs 
as well as the morphology and convert these into merger 
rates in section [6] Sections [7] and [8] are our discussion and 
conclusions. 



2 MORPHOLOGICAL PARAMETERS AND 
MERGER CRITERIA 

In this series we use the Concentration- A symmetry- 
Smoothness parameters (CAS IConselic3 20031) , com bined 
with the Gini-M2o parameters from Lotz et al] (|2004l ) and 
one addition of our own Gm. We discuss the definitions of 
these parameters in the previous papers, as well as how we 
estimate uncertainties for each. Briefly, given a set of n pixels 
in each object, iterating over pixel i with value li, position 
Xi, yi with the centre of the object at x^, yc these parameters 
are defined as: 



with r/ as the radial aperture, centered on Xc, yc containing 
pe rcentage / of the hght of the galaxy (see definitions of r/ 
in iBertin fc Arnouts|[l99^ : iHolwerdalfeoOSl ) . 



A 



Iiso\ 



(2) 



where 7igo is the pixel at position i in the galaxy's image, 
after it was rotated 180° around the centre of the galaxy. 

S„,|/(z,j)| ' 

where Is is pixel i in a smoothed image. The type of smooth- 
ing has changed over the years. We chose a fixed 5" Gaussian 
smoothing kernel for simplicity. 

The Gini coefficient is defined as: 



G = 



-E,(2i 



1)1. 



(4) 



In{n - 1) " 

where the list of n pixels was first ordered according to value 
and / is the mean pixel value in the image. 



Mao = 



log 



Mtot 



for E,7, < 0.2Jt„t 



(5) 



where Mi is the second order moment of pixel i; Mi — li x 
[{x — Xc)^ + {y — Vc)^]- Mtot is the second order moment 
summed over all pixels in the object and M20 is the relative 
contribution of the brightest 20% of the pixels in the object. 
Instead of using the intensity of pixel i, the Gini parameter 
can be defined using the second order moment: 



Gm ~ 



Mn{n - 1) 



E,(2i 



1)M., 



(6) 



These parameters trace different structural characteris- 
tics of a galaxy's image but these do n ot span an orthogona l 
parameter space (see the discussion in lScarlata et aLlbOOTl ). 
Originally, the above parameters were envisaged to classify 
the morphologies of galaxies but it was soon realized that 
a sub-space of the param eters is oc c upied by gravitation- 
ally in teracting late-types. IConselic^ (|2003l ) and iLotz et al] 
I2OO4I) introduced several different criteria for the selection 
of merging systems in ter ms of the CAS an d Gini-Mao pa- 
rameters For optical data, IConselic3 (|2003| ) define the fol- 
lowing criterion: 



A > 0.38, 



(7) 



wit h some authors re quiring A > S as well. 

iLotz et al] l|2004h added two different criteria using Gini 
and M20: 



G > -0.115 X A/20 + 0.384 
and 

G > -0.4 X A + 0.66 or A > 0.4. 



(8) 



(9) 



G = 5 log(r8o/?-2o). 



(1) 



The latter being a refinement of the Conselice et al criterion 
in equation [T] 

These criteria were developed for optical morpholo- 
gies, typically observed in restframe Johns on-B or SDSS-g. 
There fore, in the third paper in this series jHolwerda et al.| 
l2011cl) . we defined several possible criteria specifically for 
the Hi perspective using the CAS-G/M2o-GAf space of the 
WHISP survey H i map sample. We defined the Gini param- 
eter of the second order moment, Gm and a criterion for this 
parameter that selected most interacting galaxies: 



4 B.W. Holwerda 



Gm > 0.6, 



(10) 



Earlier in this series, we speculated that a combina- 
tion of Asymmetry and M20 co uld well be used to se lect 
interaction in H I morphology in l|Holwerda et al.ll2011bl '). In 
iHolwerda et al. 

I (Eoilc), we defined such a criterion as: 



A > -0.2 X M20 + 0.25. 



(11) 



Finally, we also defined one b ased on Concentra tion and 
M20, following the example of the lLotz et al] (|2004l ') criteria 
(eq.|H]and[ll): 

-5 X M20 + 3. (12) 

In l|Holwerda et al.ll2011bl ). we found that this last criterion 
selected both the correct fraction of interacting galaxies and 
that it agreed most often with the previous visual identifi- 
cations in the case of individual WHISP galaxies. We will 
now explore merger rates based on the above criteria for H I 
morphology. 



3 LIMITATIONS 

Similar to other morphological selection schemes, we note 
that our approach is most sensitive to mergers involving 
at least one gas-rich late-type galaxy for the morphological 
selection and two in the case of pair selection of mergers. 
H I observations pre-select agains early-types (see Figure [2] 
and section BTTjl and morphological disturbance is sensitive 
to unequal mass mergers (cf Lotz et al. 2010b). Therefore, 
this approach is complementary to existing morphological 
identification of mergers but dissimilar enough to warrant a 
separa te estimate of time scale s. 

In lHolwerda et al.1 (|2011dh . we compared the visibility 
time scales for the above criteria in the case of mergers of 
two equal mass spirals to those of a secularly evolving spiral. 
We find that the spiral-spiral merger is visible in morpho- 
logical criteria during two stages before the final coalescence 
into the merger remnant; once during initial approach (be- 
fore the stellar disk is disturbed) and during the second pass, 
before coalescence. The total visibility time is approximately 
a Gigayear with some variance due to observation angle, dif- 
ferent treatment of feedback from star-formation on the ISM 
in the simulation, and the relative gas fraction of the spiral 
disk. We note that the timescales for selection for merging 
and isolated (passively evolving) Hi disks become the same 
for resolutions coarser than the WHISP observations used 
here, e.g., the VLA Imaging of V irgo spirals in Atomic gas 
survey (VIVA IChung et all 120091). A limitation to the sim- 
ulations used (originally from Cox et al] l2006al lbh. is that 
they are for Milky Way-size spiral galaxies only and do not 
consider minor mergers. Thus, since our approach may be 
sensitive to some minor merger scenarios, which possibly 
have much sh orter visibility time-scal es, the inferred visibil- 
ity time from iHolwerda et al] (|2011d^ should be considered 
the upper limit for Hi mor phological selection of mergers. 

In a subsequent paper (|Holwerda et al.ll2011el ) , we show 
that the H I morphology is also sensitive to ram-pressure by a 
dense intergalactic medium but that one can select against 
ongoing or recent stripping with t he Concentration index . 
For the WHISP survey we find in IHolwerda et al] (l2011c ) 
from a WHISP sub-sample and in IHolwerda et alj l|2011dl ) 




200 300 
(km/s) 



Figure 3. The distribution of rotation velocity (vrot) over the 
complete WHISP sample. 



from simulated Hi maps, that the level of contamination 
for the above parameters varies but is acceptable for large 
volume studies. For example, these still are noisy Hi maps 
in our morphological selection (Figure [6] and the Appendix 
in the electronic version of the journal). 



4 WHISP 

The dataset here are the observations done as part of 
the Westerb ork Hi Survey of Irregular and Spiral Galax - 
ies (WHISP. Ivan der Hulst eral]|200ll : lvan der Hulstll2002l '). 
WHISP is a survey of the neutral hydrogen component in 
spiral and irregular galaxies with the Westerbork Synthesis 
Radio Telescope (WSRT). It has mapped the distribution 
and velocity structure of Hi in several hundreds of nearby 
galaxies, increasing the number of H I observations of galax- 
ies by an order of magnitude. The WHISP project provides 
a uniform database of datacubes, zeroth-order and velocity 
maps. Its focus has been on the structure of the dark matter 
halo as a function of Hubble type, the TuUy- Fisher relation 
and the dark matter content of dwarf galaxies ([Swaters et alj 
l2002al : ISwaters fc Balcells|[200^ : iNoordermeer et alj|2005al '). 
Until the large all-sky surveys with new instruments are 
completed, WHISP is the largest, publicly available data- 
set of resolved H i observations. We compiled a catalogue of 
basic data, obtained the highest available H I column density 
maps and estimated the representative volume of WHISP. 



4.1 WHISP Basic Data 

Basic data for the WHISP sample came from the UGC cata- 
logue, updated ft-om HyperLEDA (|Paturel et al.ll2003al ). We 
used updated positional data, preferring, in order, 2MASS 
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later lowered to 20 mjy. Observation times were typically 
12 hours of integration. The galaxies satisfying these selec- 
tion criteria generally have redshifts less than 20000 km/s 
{z < 0.07). 

The observational criteria (see above) are in effect a 
selection against early type galaxies (preferring spirals and 
irregulars, F igure [21), and a prefe rence for galaxies below cz— 
5000 km/s (|Noordermeeill2006l . chapter 2). Figure [2] shows 
a histogram of the Hubble types in WHISP and the same 
volume in the UGC. There is a preference for later-type 
galaxies but no exclusive selection; only a few percent of the 
early types are selected and ~ 10% of the later types. 

The WHISP data were retrieved from the "West- 
erbork On the Web" (WOW) project at ASTRON 
( http : / / www ■ astron . nl/wow/ [ ) . We use the column den- 
sity maps with the highest resolution available (~12" x 
127sin(5)). 



Figure 4. top panel: The distribution of recessional velocities 
of the WHISP galaxies and those in the part of the UGC that 
conforms to the selection criteria for WHISP. bottom panel: 
percentage of the UGC catalogue observed in WHISP as a func- 
tion of redshift. The fraction of WHISP galaxies drops off with 
redshift. The exponential fit to this drop-off is the dashed line, 
used to estimate the representative volume of the WHISP survey. 



l|Kleinmann et al.ll 19941 ), the updated Upps ala Galaxy Cat- 
alogue positional d ata jCotton et aU 19991 ') . the Principal 
Galaxy Catalogue (|Paturel et all l2003bl). the o riginal Up- 
psala Galajcy Catalogue positions (iNilsonlfigT^ ) and lastly 
the compilation of coordinates internal to HyperLEDA. The 
major and minor axis, came from the same catalogues in the 
same order. To define a sufficient sized area around the H I 
disk, we multiplied the major axis with a factor seven. This 
is to speed up computation and leave out a galaxy's com- 
panions in the column density maps. 

For the morphological information we again relied 
first on the 2MASS catalogue and secondly on the Upp- 
sala Galaxy Catalogue and lastly on any information in 
HyperLEDA. The re dshift information is primarily from 
ISprineob et all (|2005l ) for many galaxies with the remain- 
ing ones filled in from a myriad of sources in HyperLEDA. 
Figure [2] shows the distribution of Hubble types in the UGC 
and WHISP catalogue: there is a clear preference for late- 
types in WHISP. 

We also obtained HyperLEDA values for the rotational 
velocity (vrot)- Figure [3] shows the distribution of Vrot over 
the WHISP sample: WHISP selection prefers smaller («rot < 
120 km/s), and more nearby systems (Figure |4]). 



4.2 WHISP Column Density Maps 

The WHISP observation targets were se lected from t he Up- 
psala General Catalogue of Galaxies (|Nilson|[l973 ). with 
blue major axis diameters > 2^0, declination (B1950) > 
20 degrees and fiux densities at 21-cm larger than 100 mJy, 



4.3 The WHISP Volume 

A definition of the WHISP volume is not straightforward as 
WHISP was not meant as a complete volume-limited sam- 
ple of galaxies. In the case of a blind Hi survey, the esti- 
mate of the volume sampled is complicated by the detection 
function of galaxies in the observations which depends on 
the bandwidth, frequency resolution and thresh old used in 
the survey (see lZwaan et al.lll997l : IZwaanll200d . chapter 3). 
However, since WHISP is a targeted survey from an existing 
optical catalogue, we can compute the volume represented 
by the optical catalogue (UGC), estimate what fraction of 
the UGC the WHISP catalogue represents and thus what 
fraction of the UGC volume is representative of the WHISP 
survey. 

Naively, the volume covered by the UGC with a dec- 
lination over 20° and cz < 5000 km/s (r=68.5 Mpc) is: 
VuGc{S > 20°) = (27r/3) r^/i = (27r/3) x 1 - sin(5) = 
4.43 X 10^ Mpc^ Of the 8147 galaxies in the UGC, 339 galax- 
ies are in WHISP; 4.17 % of those in the volume. However, 
to equate the WHISP volume to 4.17 % of the UGC volume 
(18473 Mpc'^) would be simplistic as there is a bias towards 
nearby galaxies in the WHISP selection (see Figure |4]). 

Figure [4] shows the distribution in redshift of the 
WHISP sample and the total UGC sample (5 > 20°), as 
well as the percentage of the UGC galaxies in WHISP. We 
fit an exponential distribution to the fraction and obtain the 
radial weighting function for the volume of the UGC corre- 
sponding to the WHISP sample: 'w{r) = foe~^^'^ with fo 
= 5.7 % and h = 1947 km/s = 27 Mpc. To compute the 
WHISP volume we integrate over radius, weighting the ra- 
dius with the above function: V = J 27vrhw{r)dr = 6835 
Mpc^ 1.5 % of the UGC volume. We will use this volume 
for our computation of volume merger rates further in this 
paper. Because the WHISP survey was never meant to be a 
volume-limited estimate, this estimate of the representative 
volume should be treated with caution. Fortunately, the fu- 
ture planned Hi surveys with ASKAP and APERTIF are 
set to be volume-limited. 
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Table 2. The galaxies in WHISP with one or more compan- 
ions in the data-cube. Quah fiers of intera ction (Int?) from either 
iNoordermeer et alj ll2005bl ) (NM05) or ISwaters et"al] l l2002bl ') 
{SW02). 



Galaxy 


Companions 


Int? 


Ref 


UGC 624 


2 


y 


NM05 


UGC 1437 


1 




_ 


UGC 2141 


1 


n 


NM05 


UGC 2154 


multiple 


y 


NM05 


UGC 2459 


1 




_ 


UGC 2487 


2 


n 


NM05 


UGC 2916 


2 


y 


NM05 


UGC 2941 


1 


y 


NM05 


UGC 2942 


1 (UGC 2943) 






UGC 3205 


3 


n 


NM05 


UGC 3382 


1 


n 


NM05 


UGC 3384 


1 






UGC 3407 


3 


y 


NM05 


UGC 3426 


1 


y 


NM05 


UGC 3546 


1 


n 


NM05 


UGC 3642 


1 


y 


NM05 


UGC 3698 


1 


n 


SW02 


UGC 4458 


1 


y 


NM05 


UGC 4666 


1 


n 


NM05 


UGC 4806 


multiple 


_ 


_ 


UGC 5060 


1 


n 


NM05 


UGC 5935 


multiple 


y 


SW02 


UGC 6001 


1 


n 


NM05 


UGC 6787 


1 


y 


NM05 


UGC 7183 


1 






UGC 7353 


1 






UGC 7506 


1 


n 


NM05 


UGC 7989 


1 


y 


NM05 


UGC 8271 


3 


y 


NM05 


UGC 9642 


1 






UGC 9858 


1 






UGC 10791 


2 






UGC 11283 


1 






UGC 11951 


1 


y 


NM05 


UGC 12815 


multiple 


y 


NM05 



5 MERGER FRACTION 

There are two ways for us to estimate the merger fraction of 
the WHISP sample: by counting the number of close pairs 
or to count the number of disturbed looking galaxies. 

We should note that in the lowest mass range (M < 
10^° Mq), the observed merger fracti ons are very high 



for r e dshift range z= 



201C ; iKartaltepe ct al 
200g ; [Conselicc ct al. 2009 



0. 2-1.2 (oc 
~ 200j|: 



10% iBridge et al.1 |2007| 



Lotz e t al." '20081: iLin et al.1 
2009). Because the 



Jogee et a l. 

WHISP selection prefers nearby, irregular and smaller sys- 
tems (Figures 31 [5] and [31) , one can expect a high fraction of 
them to be merging. 



5.1 Galaxy Pairs in WHISP 

There are several galaxies that have a close companion in the 
Hi datacube. Each datacube is a single WRST pointing (10' 
X 10') with a bandwidth of 320, 680, 1280 or 2560 km/s, de- 
pending on the velocity resolution used. While this is not the 
typic al selection criterion for pair selection (see lPatton et all 
I2OOOI ). we could use it as such since pairs are selected for 



proximity on the sky and in redshift. In the full WHISP 
catalog, there are 35 galaxies with one or more companions 
in the WHISP cube. Naively this translates to a close com- 
panion and hence merger fraction {fmgr) of ~ 10 % of the 
WHISP sample. The merger fraction based on the close pairs 
depends on how many of those galaxies with companions one 
would consider merging. Typically, the velocity difference 
is taken to be less than 500 km/s to constitute a merging 
pair, so the datacube criteria are not stringent enough. If we 
go by the merger qualifiers from lSwaters et aU (|2002lj ) and 
INoordermeer et al.l (|2005lJ ). 10 of the 24 galaxies they clas- 



sify and who have companions are not merging (68% success 
rate, see table[2l). So the real merger fraction of the WHISP 
catalogue is closer to ~7 %, whi c h puts it close to the loca l 
valu es from iPatton et al.l (|l997l ). iDe Propris et al.1 (|2007al ) 
and iDarg et al. I (|2009l ) (See Figure [1]) for the local volume. 



5.2 WHISP Merger Fraction from Hi Morphology 

In iHolwerda et al.l (|201ld ). we identified the part of mor- 
phology parameter space that contains a representative 
number of the merging galaxies in a subsample of the 
WHISP database for which we had visual classifications 
of inte rac tion using the Himaps from either jSwaters et al.l 
l2002bl ) or INoordermeer et all (|2005bl ). Based on a plot sim- 
ilar to Figure [5l we concluded that criteria based on Asym- 
metry, M20, Concentration and Gm selected the correct frac- 
tion of interacting galaxies in a given sample (eg. 1101 - H2p . 
Especially in the case of the Concentration- A/20 selection, we 
obtained not only the same fraction of interaction but this 
criterion agreed with the majority of the visual classification 
of the H I map in individual ca ses. And in this paper's com- 
panion (|Holwerda et al.ir2011dl ). we explored how long both 
these criteria and those from the literature (eq. 0® select 
mergers in by their H I morphology. 

We can now apply these selection criteria to the full 
WHISP sample. The values for the morphological param- 
eters of the full WHISP sample are listed in Table ?? in 
the Appendix {electronic version only). Table O summarized 
our results for morphological selection for each of the six 
criteria: the fraction of the total WHISP sample selected, 
the resulting volume d ensity, the visibility timescale from 
iHolwerda et al.l (|2011cJ ), and the computed volume merger 
rates. For comparison, it also shows the values for the merger 
selection based on close pairs computed above. Figure [SI 
shows the parameter space highlighting those selected by 
the Concentration- M20 criterion. We excluded those galaxies 
with A=l, as this extreme value is indicative of an incorrect 
central position {xc,yc) of the WHISP galaxy from Hyper- 
LEDA. Starting with the best performing selection criterion 
(eq. [ni, we find 45 galaxies out of the 339 in the WHISP 
catalog are interacting, or 13 %. The other selection crite- 
ria select much higher fractions. The next best performing 
criterion (Gm) selected mergers very cleanly in the N-body 
simulations but its timescale appears to be very resolution 
sensitive. The Concentration- M20 criterion selected these 45 
galaxies based on their Hi morphology but we do not ex- 
pect each to be a merger individually (see Figure [6l for some 
examples from the selection. All the H icontours maps over- 
plot on 2MASS-K images are shown in Figure Al in the 
Appendix electronic version). Close to the selection criterion 
(dotted line in Figure [5l panel IX), individual measures may 
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Morphology Parameterspace 

• non-interacting 

• interacting 




0.2 0.4 0.6 0.8 1 0.5 0.6 0.7 0.8 -3 -2 -1 00.4 0.6 0.8 
A G G„ 



Figure 5. The distribution of morphological parameters, Concentration (C), Asymmetry (A), Gini (G) and the contribution to the 
second order moment of the brightest 20% of pixels (M20), and the Gini coefficient of the second order moment of the pixels (Gm)- 
Merger selection criteria from the literature are marked with da s hed lin es in panel II (equation [Sj , panel IV (equations [7] and |9ll , and V 
and VI (equation O . Our selection criteria from iHoIwerda et al.l l l201ld ) are marked with dotted lines; the Gm criterion in panels I, III, 
VI and X (equation llOp . the A-M20 criterion in panel V (equation and the C-M20 criterion in panel IX (equation [12} . Those objects 
selected by this last criterion (additionally requiring that Asymmetry is not extreme; A!=l) are marked in the plot to illustrate. WHISP 
morphological values are in Table ?? in the electronic version of the manuscript 

scatter in an d out of the selection. H owever, based on our in the earhest stages of the mergei0. From the morpholo- 
experience m iHolwerda et~aD (|201ld ). the fraction of merg- gies of the selected galaxies (Figure [6] and the Apppendix 
ers in WHISP is correct and individual galaxies are likely 
to be interacting viewed in the H I perspective. Often their 

optical appearance may be still undisturbed as these are ^ee also the Hi Rogues gallery for many examples of 

disturbed H I morphology in normal appearing spirals 
http : / /wB M . nrao . edu/astrores/HIrogues/RoguesLiving. shtml 
[Hibbard et al.l 1I2OOII I. ' 
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Figure 6. Three random examples of the galaxies selected by the Concentration- M20 criterion (eq. I12I I. The grayscale image is the 
2MASS K-band image and the contours are WHISP column density contours at 2.5, 5, 10 and 20 X 10^'^ atoms/arcsec^. . The full set 
of galaxies selected by this criterion is shown in the Appendix {electronic version only). 



in the electronic version of this paper), one can see that 
some noisy maps are still selected as well. A more uniform 
(pipeline-reduced) survey will suffer less from these spurious 
selections. 

The merger fraction we find from H I morphology selec- 
tion is highe r than other autho rs find for the local universe; 
for example, iDarg et al.l ll2009l'l find 1-3% of all galaxies in 
SDSS to be merging and lPatton et al" I (|2002f ) similarly find 
only a few percent from galaxy pairs. However, our fraction 
is similar to those found at slightly higher redshifts (z~0.1, 
see Figure[T]). As pointed out in section 4, the H I morphology 
selection is likely sensitive to some minor merger scenarios 
as well. Minor mergers are expected to dominate the num- 
ber of ongoing mergers and could in part explain our higher 
fraction. 



0.7 X 10~^ mergers Mpc~^ Gyr~^. Merger rates based on 
H I morphology can use a variety of selection criteria (Table 
[31) and the merger rate from the Concentration- M20 criterion 
is Rmgr(morph) = 6.8 x 10"^ mergers Mpc"^ Gyr"\ 

In comparison to the fraction found from galaxy pairs 
(7%) with the number found from morphology (13%), both 
seem to be in reasonable agreement with the morphologi- 
cal selection on the high side. This is in line with the dis- 
crepancy found between pair selection and morphological 
selection as can be seen in Figure [T] After factoring in the 
relat ive timescales a merger is visible as a clos e pair (0.5 
Gyr iPatton et all I2OO0I : iPatton fc Atfieldl |2008|) or above 
the Co ncentraton-M2o criterion (~1 Gyr; Holwerda et al.l 
l2011dl ). the volume merger rates agree very well: 7 and 
6.8 X 10~^ mergers Gyr~^ Mpc~^ for paired and morpho- 
logically disturbed Hi disks respectively. 



6 WHISP VOLUME MERGER RATE 

From the number of galaxies with a companion or the 
number of disturbed looking galaxies in a given volume 
(ric and ndiat respectively), one can calculate the vol- 
ume merger rate (Rmgr), provided one has an estimate 
of the merger time scale (Tmgr), the merger rate from 
pairs; Rmgr (pairs) = Uc/Tmgr (pairs) or the merger rate 
from morphology Rmgr (morph) = ndist/Tnigr(morph). From 
iHolwerda et"alT (|2011dl ) . we have an estimate of the mean 
merger time scale with some variance due to differences in 
merger conditions (type of feedback physics in interstellar 
matter, type of encounter, and gas ma sses of the d isks) a nd 
perspective (face-on versus edge-on). iLotz et aL I (|2010al lbl) 
note similarly that time scales depend on mass ratio and 
gas fraction for optical morphological selection. 

Mergers were on average visible for 40% of the 2.5 Gyr 
of the merger simulation, making our typical ti me scale 
Tmgr (morph) ~ 1 Gvr.fsee iHolwerda et al.i r2011dl . this pa- 
per's companion paper.), very similar to those used in th e 
literature for morpholog ical selection. IPatton et al.l l|200Cf ): 
iPatton fc Atfieldl (|2008l ) use a merger time scale for pairs 
of Tmgr (pairs) — 0.5 Gyr. The volume represented by the 
WHISP sample was computed above f ^43)) as 6835 Mpc^. 

Following our simple merger fraction of 7% from the 
number of WHISP galaxies with companions in the dat- 
acube, we obtain a volume merger rate of Rmgr (pairs) = 



7 DISCUSSION 

The WHISP sample represents only a very small volume of 
the Universe and the resulting merger fraction and rates are 
uncertain as a result of that. However, the reasonable con- 
sistency with muc h larger samples such as the Sloan Digital 
Sky Survey (e.g., iDarg et al.ll2009l . see Figure [T| are cause 
for optimism as Hi morphology as a tracer of the merger 
fraction and rate of galaxies. 

Volume merger rates in the literature for the local Uni- 
verse vary somewhat with sample and survey. Masie di et al.l 
(|2006l ) finds for luminous red galaxies in SDSS a volume 
merger rate oiRm ov = 0.6 x 10"* Gpc~^ Gy r~^ = 1-7 x 10"" 
h^ Mpc"^ Gvr~^ lDe Propris et aU (|2007bl ') found for galax- 
ies of all types in the Millennium Galaxy Catalogue (MGC) 
a volume m erger rate of Rmav = 5.2 ± 1.0 x 10~" h^ Mpc~'^ 
Gyr-\ and IPatton fc Atfieldl l|2008l ) find a volume merger 
rate for all galaxy types based on SDSS and MGC of Rmgv 
= 1.4 ± 0.1 x 10"'* h^ Mpc~^ Gyr~^ for major mergers. In 
contrast, we find a volume merger rate of Rmgr (morph) = 
2 X 10"^ mergers h^ Mpc"^ Gyr~^ (h = 0.73 or ^0=73 
km/s/Mpc), an order of magnitude more than those above. 
Since our merger fractions are similar within a factor two to 
those in the literature for the local Universe, the issue for 
the volume merger rate would either have to be the inferred 
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Table 3. Interaction fractions, merger visibility time and merger rate for the WHISP sample based on different morphological selection 
criteria 



Criterion 


Mergers 


fmrg 


^mrg 


T 

mrg 


^mgr 




# 


(%) 


mergers Mpc~^ 


Gyr 


mergers Gyr~^ Mpc~^ 


pairs 


15 


7 


0.0035 


0.5 


0.007 


A > 0.38 


221 


65 


0.032 


1.85 


0.017 


G > -0.115 X M20 + 0.384 


178 


53 


0.026 







G > -0.4 X A + 0.66 


235 


69 


0.034 


0.15 


0.23 


Gm > 0.6 


81 


24 


0.012 


0.80 


0.015 


A < -0.2 X M20 + 0.25 


151 


45 


0.022 


0.9 


0.025 


C> -5 X M20 + 3 


45 


13 


0.0066 


0.97 


0.0068 



WHISP volume, the timescale or a bias in the selection of 
galaxies. 

The merger time scale is unlikely to be the issue. The 
visibility time of the merger starts earlier in H I than in 
the stellar perspective but it is not substantially different 
from what other authors have found. Substituting any other 
visibility time-scale from the literature for morphological or 
pair selection would not reduce the merger rate (selection 
times are typically less or equal to ~ 1 Gyr). We are more 
sensitive to minor mergers and the implied shorter visibility 
time-scales but this is unlikely to be an order of magnitude 
effect. 

Alternatively, we may have to consider the possibil- 
ity that the morphologically disturbed galaxies are not all 
gravitationally disturbed but may suffer from effects unique 
to the gas perspective, for example ram-pressure strip- 
ping affecting the appearance of the H I disk. Compared to 
the observed fra ction of mergers from other sources (e.g., 
iDarg et al.|[2009l . for the SDSS), this is of order a factor two 
discrepancy. The agreement between volume merger rates 
from WHISP from the pairs and morphology contradict this 
however. 

The WHISP volume estimate is a likely source of the 
discrepancy as it is the most uncertain of our numbers. How- 
ever, even with the naive, larger estimate of the WHISP vol- 
ume (18 X 10^ Mpc"^), this is only a factor three, not an order 
of magnitude. 

The WHISP selection process favours late-type galaxies 
(Figure [5]) and local small irregular galaxies (to complement 
the spirals) and it was never intended as a volume-limited 
sample. Hence an intentional or unintentional selection bias 
may well have been introduced. Mergers identifiable by their 
morphology are more likely to happen to the gas-rich late- 
types and the irregulars are confined to a local -smaller- 
volume, and many of them will be tidally affected. The H I 
perspective is likely to be more sensitive to unequal mass 
mergers as these can be identified much more readily (the 
contrast in gas surface d ensities is not as gre at as it is in stel- 
lar surface brightness). iLotz et akl (|2010bl) points out how 
one expects a much higher merger fraction in lower-mass 
systems and lPatton fc Atfieldl (|2008l ') points out that a fac- 
tor two discrepancy can easily be expected if lower mass 
systems are included in eve n a p air statisti c al an alysis. In 
addition, iLotz et al.l (|2010al ') and IConselicd (|2009l ') identify 
gas-rich mergers are the most easily identified by their mor- 
phology. 

Therefore, we suspect our result points to a higher 



merger fraction and volume merger rate for spirals and ir- 
regulars in the local Universe, and less to a gross error in our 
WHISP volume, merger fraction or the merger time-scale. If 
this is the case, merger fractions and rates for at least these 
types may not evolve with redshift as dramatically as pre- 
viously thought. However, to confirm this, one would need 
a volume-limited large Hi survey with sufficient resolution 
and sensitivity for both morphological selection as well as 
accurate pair identification. 

The Hi perspective can be reliably used in the local 
Universe(z~0), where a spatial resolution can be achieved 
in large, all-sky surveys (e.g., WALLABY on ASKAP, Ko- 
ribalski et al. in prep, and WNSHIS on WSRT, Jozsa et al. 
in prep.) but any morphological identification of interact- 
ing gas disks at higher redshift will have to wait for SKA. 
The resolution of the Pathfinder instruments (MeerKAT 
and ASKAP) may well be enough to identify close pairs in 
the proposed deep H isurveys (DINGO on ASKAP, iMeveJ 
(2009), Meyer ct_al. vn^prep, and LADUMA w ith M eerKAT, 
iHolwerda fc BlvthI (|20ld '): iHolwerda et al.l (|2011a|), Holw- 
erda et al. in prep.). 



8 CONCLUSIONS 

In this paper, we explored the merger fraction and rate based 
on the Hi observations of the WHISP sample of galaxies. 
The sample is still small compared to other local references 
based on f or instance the SDSS o r the Millennium Galaxy 
Catalogue (|De Propris et al]|2005l ) but provides us with an 
indication how well the Hi surveys of the near future will 
perform in this respect. From the quantified morphologies 
of the WHISP column density maps , we conclude: 

1. The merger fraction in the WHISP sample is 7% based 
on pairs, and 13 % based on disturbed morphology. These 
percentages are consistent if one takes into account how long 
a merger is visible as a close pair of galaxies and how long 
as a morphologically disturbed H I disk. 

2. Assuming the representative volume of the WHISP 
sample is 6835 Mpc^, and a merger visibility time scale 
of 1 Gyr, the merger rate for our selection criterion is 
Rmgr(morph) — 6.8x 10"'^ mergers Gyr~^ Mpc~^ in the local 
Universe, very close to the value of Rmgr(pairs) = 7 x 10"'^ 
mergers Gyr~^ Mpc~'^ for galaxy pairs in WHISP. 

3. While the WHISP merger fractions and especially rates 
mutually agree, the merger rates are much higher than those 
reported in the literature. Selection effects in the WHISP 
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survey, preferring dwarf and irregulars, rather than a gross 
error in the WHISP volume, could well account for the dif- 
ference as well as variation in the quality of the WHISP 
maps across the sample. Upcoming, volume- limited H I sur- 
veys should provide an accurate measurement of the local 
merger rate from both Hi morphology as well as close pairs. 



9 FUTURE WORK 

The 21 cm window on the Universe is set to revolutionise 
our understanding of the merger rate of spiral and irreg- 
ular galaxies as three independent measures of merging or 
gravitational interaction are available in the data: Hi mor- 
phology, kinematic signatures in Hi of interaction (e.g., lop- 
sidedness of the profile, non-circular motions and an irregu- 
lar velocity field) and the easy detection of physically close 
companions in the Hi data-cube. 

Planned Hi surveys with the SKA Pathfinder instru- 
ments, include an all-sky survey (WALLABY with ASKAP) 
and (WNSHIS with WSRT/APERTIF), and a medium deep 
survey to z ~ 0.4 (DINGO) with ASKAP, and an extremely 
deep survey (LADUMA) with MeerKAT. Combined, these 
will revolutionise the volume probed with H I and help shed 
light on the merger fraction using all three tracers; morpho- 
logical, dynamic and pair identification in the local Universe 
surveys, and dynamical and close pairs identification our to 
higher redshift. 
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